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We calculate the relativistic constraint equation which relates the curvature perturbation to the 
matter density contrast at second order in cosmological perturbation theory. This relativistic "sec- 
ond order Poisson equation" is presented in a gauge where the hydrodynamical inhomogeneities 
coincide with their Newtonian counterparts exactly for a perfect fluid with constant equation of 
state. We use this constraint to introduce primordial non-Gaussianity in the density contrast in the 
framework of General Relativity. We then derive expressions that can be used as the initial condi- 
tions of N-body codes for structure formation which probe the observable signature of primordial 
non-Gaussianity in the statistics of the evolved matter density field. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Our knowledge of the statistics of the galaxy distribution relies upon the vast amount of data obtained by increas- 
ingly large galaxy surveys [HI]. Among other goals, analysis of the galaxy field allows us to indirectly probe the 
distribution of the underlying dark matter on non-linear scales (e.g.,Refs. On the theoretical side, in order to 

understand the physics that governs the observed galaxy field, large numerical codes are developed to simulate the 
evolution of matter inhomogeneities that have formed large scale structure (LSS). This huge task is usually split in 
two stages. In a first stage, semi-analytical methods are employed to account for the early evolution of fluctuations 
in the weakly non-linear regime. At the same time, the inhomogeneous in the continuum matter field are related to 
a discrete distribution of point masses, thus implementing initial conditions for numerical codes. In a second stage, 
typically at redshifts z ~ 50, N-body codes evolve inhomogeneities in the strongly non-linear regime up to the present 
day. As Newtonian N-body codes continue to improve in resolution and volume (e.g., Refs. [ZriS]), the implementation 
of realistic and accurate initial conditions is increasingly important. 

Historically, the initial conditions for N-body simulations have been generated by using the Zel'dovich approximation 
|10j . which establishes the correspondence between the matter density fluctuation of standard perturbation theory, 
and the displacement of mass particles in a grid. Despite its linear nature, this represents an improvement over 
standard perturbation theory, since it takes advantage of working in Lagrangian coordinates [5J 111] . The caveat to 
this approximation is that it accounts only for the early non-linear evolution of density fluctuations, and in particular, 
it employs a linear Poisson constraint, which is used to express the density contrast, <5n, in terms of the gravitational 
potential, </>n, that is 

V 2 (/> N = 4vrGp a 2 '5N. (1.1) 

An improvement to this approximation is second-order Lagrangian perturbation theory (2LPT), which generates 
initial conditions taking into account non-linearities in Lagrangian coordinates. This has been shown to be more 
precise and avoids transients present in the Zel'dovich approximation |12l 113) . Since 2LPT takes into account non- 
linearities, the fact that the gravitational instability is non-local is manifest in corrections to Eq. ( |1.1[ ) given by tidal 
effects at non-linear order |14j . With the matter density fluctuations at non- linear order under control, recent studies 
have used 2LPT to include primordial non-Gaussianity in the matter fluctuations |15H17j . 
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These and other semi-analytical approximations to the early evolution of inhomogeneities, however, rely on New- 
tonian physics, thereby ignoring the effects of General Relativity (GR). Cosmological inhomogeneities are well de- 
scribed by Newtonian dynamics only when the modes of the perturbations lie well inside the horizon, i.e. when their 
wavenumber is k 3> %, with % denoting the Hubble parameter in conformal time. Yet, the initial conditions for these 
approximations come from much earlier times -typically the epoch of decoupling-, when some of the scales of interest 
are comparable to, or even larger than, the cosmological horizon. Therefore, relativistic effects are important and 
should be taken into account when setting the initial conditions to simulations of structure formation. 

Recent studies demonstrate the importance of GR in the analysis of large scale structure. Some have contrasted 
relativistic and Newtonian fluctuations by the identification of dynamical equations [IE1 E5]. This provides cor- 
respondences between Newtonian fluctuations and relativistic perturbations in a specific gauge at linear order in 
perturbation theory. Additionally, Ref. [20] extends this correspondences to second-order perturbations. In this 
way, the equivalence of the dynamical equations is established for the restricted case of pressureless matter and 
neglecting the decaying mode of perturbations. A major motivation to study this correspondence is to discriminate 
primordial non-Gaussian fluctuations from non-Gaussianitics induced by the non-linear dynamics of GR. In search 
of observational signatures, Ref. |H] studied the effects of relativistic non-linear fluctuations in the halo bias and 
subsequently its signature in the spectrum of the galaxy distribution (see also [22] ). 



In this paper we present the Poisson equation at second order in the framework of relativistic cosmological pertur- 
bation theory |23H27| . Previous studies have explored this constraint for the limit of a dust universe at small scales 
[2"0] (in this case the linear equation is recovered), and for a ACDM universe at large scales [5S]. Instead, our 
analysis yields the Poisson constraint equation in terms of relativistic perturbations that find a direct correspondence 
with Newtonian inhomogeneities, and without approximations. Furthermore, we extend the constraint to the case 
of a general perfect fluid. As an example, we subsequently use our result to express the primordial non-Gaussianity 
in terms of the dark matter density field in equations that include all the relativistic effects. We present results in 
the form of kernels for the non-linear variables, a form customarily used in the formulation of initial conditions of 
numerical simulations. 

The paper is organised as follows. In the next section we explicitly show how to construct the Poisson equation 
from Einstein's field equations combining variables in two gauges for linear perturbations. In Section III we repeat 
the procedure for the second order variables and arrive at a GR version of the Poisson constraint valid for any perfect 
fluid including entropy (or non-adiabatic pressure) perturbations. In Section IV we apply the constraint to the case 
of matter perturbations in a flat universe dominated by pressureless matter and show how to include the primordial 
non-Gaussian corrections in the Poisson equation. We conclude in Section [V] discussing the relevance of our result to 
the initial conditions of numerical simulations. 



II. THE POISSON EQUATION AT FIRST ORDER 



A. Background and first-order equations 



In cosmological perturbation theory, considering scalar perturbations of the metric yields the following line element, 



J 2 2/ 

as = a 



(rj) [ - (1 + 2<j))drj 2 +2B ji dx i dr)+ (1-2^ + % diW}, (2.1) 



where <fi is the lapse function, ip is the curvature perturbation, and B and E make up the scalar shear. All these 
quantities are function of Cartesian coordinates, Xj, and conformal time, rj. Perturbations are then expanded order- 
by-order in a series as, e.g., <fi = <j>x + g<^2 + • ■ • • In order to define the expansion uniquely, and as an excellent 
approximation to observations, the first order quantities are chosen to have Gaussian statistics. 

In the background the metric represents the Friedmann-Lemaitre-Robertson- Walker spacetime. The homogeneous 
equations are the familiar Friedmann and continuity equations: 

m 2 = 8irG Po , (2.2) 

p' Q = -m{p + Po ), (2.3) 

where the prime denotes a derivative with respect to conformal time and a subscript zero denotes the background, 
homogeneous quantities. 

The fluid equations are derived from the vanishing covariant derivative of the energy momentum tensor. At first- 
order in perturbation theory, the energy conservation dictates the evolution of the density perturbation dpi, 

5p[ + m (5 Pl + 5Px) = {p + P ) [3V4 - V 2 (E[ + v x )] , (2.4) 
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where the scalar velocity potential is v(x J 1 t) and the energy density and pressure fluctuations are denoted by Sp and 
6P, respectively. We define the Laplacian as V 2 = djd^ . Note that no gauge has been specified here. In order to 
obtain the corresponding equation for the evolution of the velocity, we define V = B + v, and write the momentum 
conservation equation, which at first order is 



Vl l +H(l + c 2 s )V hl + 



5Pi 



Po+Po 



= 0, (2.5) 



where we have neglected anisotropic stresses and defined the adiabatic sound speed as c 2 = Pq/p'q- 

The Einstein field equations yield two constraint equations that are combined to derive the Poisson equation. The 
(0, 0) component of these equations yields the energy constraint equation 

3H(ip[ +U4>i)- V 2 (Vi + U{E[ - Bij\ = -Ai:Ga 2 5pi. (2.6) 

The momentum constraint is derived from the (0, i) component: 

ij/ lti + U4>x,i = -AixGa 2 {p Q + P )V lti . (2.7) 

This is the complete set of equations at first order without the gauge specified. The remaining Einstein equations at 
this order are related to the ones above through the Bianchi identities. 

B. Constraint in the longitudinal gauge 

In order to overcome the ambiguity in the coordinate freedom, we must specify the gauge in the above equations. 
We work in the longitudinal or Newtonian gauge [23] to recover the exact Newtonian equations. This is a shear-free 
gauge, specified by setting E( = Bg = 0. The absence of anisotropic stresses also guarantees that, in this gauge, 



ipit — $it an d Eq. (2.6| becomes 

- 3H(<p' u + H(j> u ) + V 2 -/»i£ = 4irGa 2 6 Pu - (2. 



Then, by integrating the overall gradient of the momentum constraint (2.7), we have 

<t>'u + H<t>u = -47rGa 2 (p + P )vu, (2-9) 
and combining both equations we find the first-order constraint: 

V 2 4> u = 4TrGa 2 (Sp u + p' Q v u )- (2.10) 

C. The Newtonian expression 



The quantity in parenthesis in the linear Poisson equation (2.10), is equivalent to the density contrast in two other 
gauges as we will now show. The transformation between two coordinate systems is carried through the four-vector 
Cifi = ( a iiPi,i)i so that, for example, the density perturbation at linear order is transformed as 

Jp 1 =5p 1 + ap' Q . (2.11) 

The total matter gauge is defined by a vanishing total momentum at all orders, i.e. 

VW = 0. (2.12) 

The transformation rule for V tells us that 

V 1 = v l + B 1 =v 1 +B 1 -a 1 , (2.13) 
so that in the case of a transformation from the longitudinal to the total matter gauge we have 

Qtltom = vu ■ (2-14) 
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In consequence, the density fluctuation is obtained as 

Spltoia = Spu + PqVU ■ (2-15) 

It is now straightforward to recover the Newtonian form of the Poisson equation writing 

V 2 fa e = 4nGa 2 p 5 Uom , (2.16) 

where the density contrast is, at first order, <5itom = ^/Oitom/pcb and at second order <52tom = Sp2tom/po- 

An equivalent transformation can be performed to arrive at the comoving gauge, where the shear is vanishing and 
B CO m — E com = v com = 0. In this case, one finds that a com = vg 1 just as in the total matter gauge. Thus, the matter 



density at linear order 8p\ com reproduces the expression in (2.14). The corresponding Poisson equation 



V 2 (f>u = 4irGa 2 p S lcom , (2.17) 

has been recovered in previous works [5J 129) . It has further been shown that with the same combination of variables 
(namely Si com , <f>u, vg) one can reproduce the equations used in Newtonian hydrodynamics at linear order [T51 fTO] . 
with the exception of fluids with non- vanishing pressure, and which allow for entropy perturbations |30| . At second 
order, however, the gauge transformation &%(_ —¥ <$2com include time integrals which may introduce non-local terms. 
We avoid such complication here by working with the matter fluctuation in the total matter gauge. For the sake of 
completeness, we construct the scalar /3tomi which is determined by the transformation Ei iom = En + pitom = 0, that 
is 

Atom = 0. (2.18) 

III. THE CONSTRAINT AT SECOND ORDER 

In the previous section we have shown how the linear energy and momentum constraint equations can be combined 
to obtain a Poisson equation at first order. The same procedure can be followed to write a Poisson-like constraint at 
second order, although the manipulation of terms is obviously more complicated. 

A. Second-order equations 

The energy constraint at second order in a non-specific gauge form is [31] . 

m{i> 2 ' + n<\> 2 ) + v 2 (h(b 2 - e 2 ') - ife) + w 2 b 1 (v 2 (s/ - - 2vv) 

+ B lA (H{mB^ - 2V 2 £i/ - 2(^i + <M/) - *ti?) + ZEi'Hih. - 2%Bi),« 
+ 4H(ih ~ <h.)(tyi' - V 2 (S/ - + E 1 '^{mE 1 + X -E^ - B^ 
+ V-i'(2V 2 (^i' - 2HE 1 ) - 3VV)) + fa. i (2W 2 E 1 - 3^i),< + 2V 2 fa{V 2 E 1 - 4^i) 
- l2H 2 fa 2 + (/;, ,,,/»':." + V 2 £?ijV 2 £i, J - E liijk E 1 y k - V 2 E 1 'V 2 E 1 '} 

= -AnGa 2 (2(p + P )V 1 k v hk + Sp 2 ^ , (3.1) 

while the momentum constraint is 

W 4 + Ufa , - E^M, + fa + V 2 E x )i + B 1Aj (2HB 1 + fa)? 



^(V 2 ^ - 4V>i)J - fa ti [m<h. + *l>i + * 7 {Ei - B l 

, ' + ^\ jE i,ii + Kik E hi jk ~ V-i^V 2 ^ + Hi) - V 2 faB, 



-AirGa 2 



(po + Po)(v 2<i - 2<j>x{V 1 + Bi),, - 4(^M,i - E liik vx, k )) 
+ 2{Sp l +SP l )Vi, i ]. (3.2) 
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Both equations are simplified by choosing a gauge. In the longitudinal gauge (incidentally called Poisson gauge at 
second order when one ignores the vector and tensor modes), Eq. (3.1 1 takes the form 



= -AirGa 2 po{Sp 2e + 2(1 + w)v Ut k v u ,kj , 



while Eq. (3.2) is reduced to 



(3.3) 



(ip' 2e + H(j)2t) t i + 2 i^u^u)' - %H4> u .i4>u - 24>u,i<l>u 

= -47rGa 2 p {(l + w) [v 2 e,i - 6vu,i4>u] + 2(1 + c 2 )v 1Li 8u + 2 — 6P nad iVu,i\ 
L Po > 



(3.4) 



Here w — Pg/po is the equation of state of the fluid, and the non-adiabatic pressure perturbation, <5P na di, is defined 
as 

^Pnadl = <5Pi - c 2 J Pl . (3.5) 



Following the steps of the procedure at first order, we take the spatial divergence of Eq. (3.4) and integrate with 
the inverse Laplacian operator V -2 . We obtain 



V4f + Hfoe + (tit)' - m<? u - 2V- 2 {(jy'u^u) J 



= - 4irGc 



a 2 P J{\ + w) \v 2 i - 6V~ 2 (vuj&u) j ] + 2(1 + c s 2 )V 2 (v ltJ S u ) 3 + 2— (SP nadl v U j) 4 

l L J po 'J 



We can now substitute this into Eq. (3.3) to arrive at 
,2 3 



V=Vm + 3 (<p'uY + ^V\<t>\t) + S<j>uV 2 <t>u + m U 2 U )' - 6HV 2 r - ' 



-AttG 



a 2 p {s 2 e - 3K(1 + w)v 2( + 2(1 + w)v Ul jV U , j + 6nV~ 



V,j ( 3(1 + W)4>U - (1 + C 2 )&U ~ —SPnadi 



1 

Po 



(3.6) 
(3.7) 

)].'} 



B. The Poisson equation at second order 



To write the second-order equivalent of the Poisson equation in (2.16), we must transform the density contrast to 
the total matter gauge. The transformation rule at second order is 26\ 

+ "item (Po^ltom + Po a 'ltom + ^Pu) ■ (3.8) 

The second-order a^tom evaluated in the longitudinal gauge is found with the aid of expressions in Ref. |31j . 

1 



o^tom = v 2 £ - Uv l u + ^(vfe)' - 4V 2 [v u ,j (4>u - v'u)] ! 



(3.9) 



With the aid of the background equations and the expressions for aitom in Eq. (2.14) and S Uom in Eq. ( 2.15[ ) we 
obtain the gauge transformation, 



S 2 tom = 5 2 g - 3"H(1 + w)v 2e - 3H(1 + w)vuSuom + 2vu5' lt 



+ 12K(1 + w)V~ 2 [v u ,j {<t>u + v' u )} J + 



SHw' - -H 2 {l + w) (5 + 9w) 



vie- 



(3.10) 



We substitute the Se factors at both orders into Eq. (3.7) for the total matter gauge equivalents. The final expression 
in terms of 5 tom , 4>u and vu is then 

vV 2 , + 3 {4>' u f + ^v 2 (0 2 £ ) + h4>uV 2 4> u + in {tf u )' - mv- 2 [M Uj ] _ j 



= AnGa 2 pol <5 2 tom + 6H(1 + w)v u 8u om - 2v u 5' ltoni + 2(1 + w)v u ,jVu, J + ijH 2 {l + w) {3w - 1) v\ t 



6H(l + w)V- 



vu,j ( <t>u - 2v[ 



1 + w 



1 6P 



nadl 



1 + w po 



1 J 



(3-11) 



() 



This rather long equation fulfils our first goal, to provide a Poisson equation at second order using the same variables 
employed in the structure formation studies at the Newtonian limit. It is already clear that adopting the expression 
V 2 i/'2 = 47rGa 2 p ^2 leaves out most of the ter ms of the actual second order Poisson constraint. 

To conclude this section, let us rewrite Eq. (3.11 1 in terms of the potential fag instead of fai. This will come handy 
in the next section since primordial non-Gaussianity is conventionally formulated in terms of this variable. We use 
the traceless ij component of the field equations as 



V 4 C0 2 £ - fat) = ~ 4V 4 (</> 2 £ ) + 2fa e /4> u / + 6V 2 fa e V 2 fa e t- oq>u,j 



V 2 fa! j 



+4irGa 2 Po {2(1 + w) [V 2 (v u ,jV U , j ) + iV 2 {v u V 2 v u ) + ZV 2 (v u )V 2 {v u ) - 3v u V 4 v u ] } • 



Upon substitution of this in the constraint equation, Eq. (3.11 ), and with some algebra we arrive at 



V fan 



2V 4 (^ £ ) + 7V 2 {fa e V 2 fa e ) + 3(V 2 <M 2 " ZfatfUu + 3V 2 (0 + +3«VV W - QU {^u^fai) , 

3. 



= 4irGa 2 p Q { V 2 <5 2tom + 6W(1 + w)V 2 (v u 6 ltom ) - 2V 2 (v u S' Uom ) + ^U 2 {1 + w)(3w - l)V 2 v 



2„,2 

U 



(3.12) 



+ 6(1 + to) 



v u \ 7 v u - V z {v u V z v u ) - (y*vu) + H v u ,j 



>u 



2v' u - 



1 SR 



nadl 



•"ltom 



l + W po 



(3.13) 



This constraint is valid for any perfect fluid. In the following section we show how to insert this constraint in the 
initial conditions of numerical simulations of structure formation. 



IV. NON-GAUSSIAN INITIAL CONDITIONS FOR NUMERICAL SIMULATIONS 



The Newtonian Poisson equation is used at all orders as a constraint to the initial conditions in numerical simu- 
lations. However, the above constraint is the one that provides consistency with General Relativity. Imposed at an 
initial time, this constraint is met at all times if the perturbations are evolved in the context of GR. It is therefore 
useful to write the expression we have derived in terms of variables employed in numerical simulations, namely <5itom 
and vu, evaluated at some initial time. Here we derive such an expression with the aid of the first order equations, 
Eqs. (2.9|, (2.161, and the continuity equation from Ref. [30] in terms of the chosen gauge. These help us to replace 



the time derivatives in the constraint equations. After some more algebra we obtain 



»2i 



2V 4 



'it 



) + 7V 2 {4>uV 2 4>u) + 3(V<M 2 - 3&/VVm 



AnGa 2 p { V 2 <5 2t om + 6(1 + to) 



viMvu ~ (V 2 w„) 2 - h/ 2 {v u V 2 v u ) + + w)U 2 V 2 {v 2 u ) 



+ mV 2 (v u 5 u ) + 6(1 + w)U 



24, 



c 2 - 1 

1 + W ' 



1 SPrnidl 
l+W po 



(4.1) 



We emphasise that the Newtonian counterpart of this constraint is a linear equation which includes only the first 
term at each side of the equality. All the other terms bring relativistic contributions to the Poisson equation. This 
expression can be used in the numerical simulations that set initial conditions for perturbations of in any perfect fluid 
and allowing for entropy perturbations. 

To reduce Eq. (4.1) further, we can either eliminate the density contrast <5itom or the potential fan via the first- 
order Poisson equation ( 2.16| ). This proves useful when we want to make contact with formulations like the so-called 
renormalised perturbation theory (RPT)[32 , where the initial conditions are set, order by order in Fourier space, 
via recursive relations in powers of <5in (see, e.g., Ref. [6]). To reduce things further, let us focus on the case of an 
Einstein-d e Si tter universe, a flat space-time filled by dust, i.e., where w = as well as SPi = and A = 0. In this 
case, Eq. (4.1) is reduced to 



V 4 c 



> 2f 



2V 4 (0 2 £ ) + 7V 2 fa e V 2 fa t ~ + 3(V 2 </> W ) 2 



4^Gp {v 2 ( 5 2 tom + \n 2 V 2 v 2 + 2 [3w u V 4 wu - 3(V 2 v u ) 2 - 2V 2 {v u V 2 v u ) 



6HV 2 {vu6 u ) + 6H [vu,j {2<j> 



->ltom 



)]/'}• 



(4.2) 
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With the aim of incorporating our result as an initial constraint in the formulation of non-linear initial conditions 
for numerical simulations, we transform Eq. (4.2) to the Fourier space. Additionally, as is customary in structure 
formation studies, we work exclusively with the growing mode of perturbations, where <j) U — const. It is then possible 
to write all of first order variables in terms of (5it om (as in the standard perturbation theory, c.f. Ref. [5]) with the aid 



of the first order Poisson equation (2.161 and the continuity equation at first order. Explicitly, in Fourier space, 



3 Ti 2 



T-L 

vu(k) = ^S ltom (k), 



(4.3) 



The reduced Poisson equation is 

k 4 foi(k) + ^k 2 n 2 5 2tom (k) = ^J d 3 pd 3 qS 3 D (p + q-k) 



T-L 

p 2 q 2 



X j3|p + q\ 4 + ^(p 4 + q 4 ) ~ f |P + <?| V + q 2 ) f P 2 q 2 + \n 2 \p + q\ 2 }<Wp)<W<?) , (4.4) 

where 8o{k) is the Dirac delta function. This equation represents a concrete constraint for initial conditions of 
numerical simulations, consistent with GR, and written in terms of relativistic equivalents to the gravitational potential 
and the matter density perturbation. Note that, while the linear equation of Ref. [2Qj is valid for these second order 
variables at small scales, the relativistic corrections obtained here become increasingly important as the perturbation 
modes approach the horizon scale. 

Since this constraint already carries couplings between different perturbation modes, there will be some intrinsic 
non-Gaussianity induced by this second-order correspondence. This is a known effect of GR which has recently been 
explored in the CMB through the use of second order Boltzmann codes [3"3Tf3"5] . and in the matter density field |28j . 
Here we disentangle the effect of the initial constraint from the influence of the non-linear evolution of perturbations. 
To observe the type of non-Gaussianity induced by the GR constraint, we introduce three templates that constitute 
a basis for the non-Gaussian 4>2t- These templates are also a basis to represent the initial conditions of the density 
contrast with primordial non-Gaussianity, and consistent with GR, for a given model of structure formation. 

Following the convention of |36j for the non-Gaussianity in the lapse function, the local template is, 



4>\° c = + [tic ~ (tfa)] 



(4.5) 



This is preserved in super-horizon scales since <pi is constant when the universe is filled with dust. Therefore, we can 
directly substitute the primordial <p2t in the Poisson constraint ( 4.4 ) . 



In Fourier space, the local bispectrum yields the correspondence 



< loc (h\ — - f loc 



/ d 3 pd 3 qS 3 D (p + q-k) (JL 



^ltom b)<5ltom(<?) , 



(4.6) 



and we can generate a kernel for 5 l £ 



loc 
om 5 



^ 2 ^2tom(^) 



3 1 



d 3 pd q 



1 



2 2 "D 

p z q z 



/nl C ) \P- 



5 3 D (p + q~k) 

15, 4 4, 

-tvp +q )- 



35 



\P- 



\\P 2 



q 2 ) 



15 



p¥ + ^i 2 \p + gf }S ltom (p)S ltom (q) . (4.7) 



Repeating the procedure for the equilateral and orthogonal configurations, we can provide initial conditions for c>2tom 
in a complete basis for primordial non-Gaussian perturbations. We borrow the templates implemented in Ref. [17] . 
For the equilateral configuration, this template is 



\k 4 ^\{k) = \fZ J d 3 pd 3 qS 3 D (p + q-k) 



T-L 4 

) 5ltom(p)£ltom(g) 

p z q^ 



\(P 2 + q 2 )\p + q\ 2 + i{p 3 + q J )\P + qIT'Ip + q\' 



(4.8) 



while in the orthogonal case 



\k 4 4>° 2 \\k) = I d 3 pd 3 qS 3 D (p + q-k) , 



(p)£itom(g) 



> -21\p + q\ 4 + 6(p 2 + q 2 )\p + q\ 2 + ^(p> + q^)\p + q\j x \p + q\< 



(4.9) 



Finally the complementary equilateral and orthogonal kernels for <52tom are 
|tW = / d 3 pd 3 q ^6%(p + q -k) 

x is (1 + 2fZ) \p + l\ 4 + \ (6/n1 - 35) \p + q\ 2 ( P 2 + q 2 ) 3/$> + ffl V + Q j )\P + Qlj 1 

is it; q ~) 

+ T (/ + 9 4 ) - Y p 2 q 2 + ^n 2 \p + q\ 2 U t am(p)*ltom(«), (4-10) 

|k°t r o m (fc) = / d 3 pd 3 q-^5 3 D (p + q-k) 

x 1 3 (1 + 7/S2) |p + <?| 4 - | + 35) |p + g| V + q 2 ) ^f&\p + g| V + <f)\p + qlj 1 

+ ^(/ + 9 4 ) - Y p2q2+ 2 n2 \P + l\ 2 ><Wp)<W<Z)- (4-11) 

The set of these three kernels are a complete basis for the primordial bispectrum. The kernels above show that 
the Poisson constraint yields different contributions for the local, equilateral and orthogonal configurations. The 
relativistic initial conditions can mimic non-Gaussian contributions, which peak in the local configuration. This 
relativistic intrinsic non-Gaussian imprint in the matter fluctuation has a maximum in the local template where 
/jvz, = — 1 (f l NL i s particularly relevant in studies of LSS since it is the dominant configuration contributing to the 
dark matter halo bias [37]). In the equilateral configuration, we find that the intrinsic GR imprints are = 1/2. 

The result obtained for the intrinsic non-Gaussianity is close to that in Ref. |28j which indicates a value in the local 
configuration of /j?? = — 5/3. The difference is most likely due to a choice of the variable in which the result is 
expressed. A detailed analysis of the modification of /nl separating initial constraints from the non-linear evolution 
of 82 is the subject of a forthcoming paper [38] , 

V. DISCUSSION 

In this paper we have derived the relationship between the energy density fluctuation and the curvature perturbation 
at second order in the context of cosmological perturbation theory. This Poisson equation at second order, presented 



in Eq. (4.1) in full generality for a single fluid including entropy perturbations, is expressed in terms of variables 
equivalent to an Eulerian set in Newtonian hydrodynamics. We found that the Poisson equation takes on a particular 
simple form, if the matter density fluctuation is expressed in the total matter gauge, and not the comoving orthogonal 
gauge which is used at first order. The difference of the two gauges only becomes apparent at second order in 
perturbation theory. 

As an example, we calculate the second order Poisson equation in the case of an Einstein-de Sitter universe, 



and present the result in Eq. (4.4) in Fourier space. We show how to incorporate primordial non-Gaussianity into 
the matter perturbation at second order in an equation consistent with GR. In this way, we can also quantify the 
non-Gaussianity intrinsic to GR contributions. Our results generalise the non-Gaussian kernels presented in terms 
of Newtonian physics in Ref. [17] to include relativistic terms. We show that the non-linearity of GR induces a 
non-Gaussian signature on top of the primordial value. In particular we find, in the local configuration, a value 
Inl = — 1- This is close to that obtained in Ref. [3S], with the difference likely due to different gauges being used. 

The general result, the Poisson equation at second order, and the example presented in this paper provide fairly 
simple equations that can be directly incorporated into generators of initial conditions for numerical simulations. The 
initial conditions generated in this way account for general relativistic effects in N-body codes and other numerical 
simulations of structure formation. 
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